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ABSTRACT 

We investigate the emission-line properties of galaxies with red rest-frame colors (compared to the g — r color 
bimodality) using spectra from Data Release 4 of the Sloan Digital Sky Survey (SDSS). Emission lines are 
detected in more than half of the red galaxies. We focus on the relationship between two emission lines 
commonly used as star formation rate indicators: Ha and [O II] A3727. Since [O II] is the principal proxy 
for Ha at z ~ 1, the correlation between them is critical for comparison between low-z and high-z galaxy 
surveys. We find a strong bimodality in [O II]/Ha ratio in the SDSS sample, which closely corresponds to the 
bimodality in rest-frame color. Based on standard line ratio diagnostics, most (nearly all of the) line-emitting 
red galaxies have line ratios typical of various types of Active Galactic Nuclei (AGN) — most commonly 
“low-ionization nuclear emission-line regions” (LINERs), a small fraction of “transition objects” (TOs) and, 
more rarely, Seyferts. Only ~ 6% of red galaxies have line ratios resembling star-forming galaxies. A straight 
line in the [O II]-Ha EW diagram separates LINER-like galaxies from other categories, provides an effective 
classification tool complementary to standard line ratio diagnostics. Quiescent galaxies with no detectable 
emission lines and those galaxies with LINER-like line ratios combine to form a single, tight red sequence 
in color-magnitude-concentration space. Other than modest differences in the luminosity range they span, 
these two classes are only distinguished from each other by line strength. We also find that [O II] equivalent 
widths in LINER- and AGN-like galaxies can be as large as that in star-forming galaxies. Thus, unless objects 
with AGN/LINER-like line ratios are excluded, [O II] emission cannot be used directly as a proxy for star 
formation rate; this is a particular issue for red galaxies. Lack of [O II] emission is generally used to indicate 
lack of star formation when post-starburst galaxies are selected at high redshift. Our results imply, however, 
that these samples have been cut on AGN properties as well as star formation, and therefore may provide 
seriously incomplete sets of post-starburst galaxies. Furthermore, post-starburst galaxies identifed in SDSS by 
requiring minimal Ha equivalent width generally exhibit weak but nonzero line emission with ratios typical of 
AGNs; few of them show residual star formation. This suggests that most post-starburst galaxies may harbor 
AGNs/LINERs. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: fundamental parameters — 
galaxies: ISM — galaxies: statistics 


1. INTRODUCTION 


[O II] A3727 emission line is a widely used, empiri¬ 
cal star formation rate (SFR) indicator, especially at 
high redshift when Ha moves out of the optical window 


(e.g. 

Gallagher et al.lll989HKennicuttJ|1992HCowie et al. 

199G 

19971 lEllislll997l IHammer et al.lll997l IHogg et al. 

1998 

Rosa-Gonzalez et al.ll2002l iHiopelein et al.ll2003ll. 


Although its luminosity is not directly coupled to the 
ionizing flux, and it is very sensitive to reddening and 
metallicity effects, it still can be empirically or theo¬ 
retically calibrated through comparison wit h Ha. A 
varie t y of calibrations have been developed llKen niaittI 
19921 ll998t lKewlev et al.112004 iMoiihcine et al.1120051 
Moiistakas et al.ll200,^ . enabling SFR estimations from 
[O II] to a reasonable accuracy for star-forming galaxies. 

However, many red, elliptical galaxies at z ~ 0 also 
have si g^fica^^ |0Jll^and^ o^er^line em^ sion in their 

also indicate star formation in these galaxies? It has long 
been realized that star formation is not the only possi¬ 
ble source of [O II] emission in galaxies. Active Galac¬ 
tic Nuclei (AGN - especially LINERs), fast shockwaves. 


post-AGB stars, and cooling flows might also produce 
[O II] emission. Thus, before we use [O II] as a univer¬ 
sal star formation indicator, we need to know how often 
and to what degree [O II] emission is contaminated by 
sources other than star formation, especially in red, el¬ 
liptical galaxies. 

The assumption that [O II] measures star formation 
has been fundamental to many studies. For instance, 
the lack of [O II] is often used as a criterion for post- 
starburst galaxy identification. Post-starburst galaxies, 
as their name suggests, are galaxies that underwent a 
strong recent star formation epoch but have stopped 
forming stars. Their spectra can be modeled by a com¬ 
bination of an old stellar population (similar to that of a 
K giant star or an early-type galaxy) a nd a young popu¬ 
lation which is dominated by A stars (|Dress]er fc Giimil 
[iM3). Therefore, these galaxies are commonly known 
as ‘K-l-A’ or ‘E-l-A’ galaxies. With time, these galaxies 
will display an early-type galaxy spectrum after all the A 
stars die out in 1 Gyr. Thus, the study of these galaxies 
is important for understanding galaxy evolution, given 
the current uncertainty in how elliptical galaxies form. 
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The identification of these post-starburst galaxies 
requires a total lack of emission lines, to be cer¬ 
tain that star formation has truly stopped. Most 
works on post-starbursts le.Ef.. IPressler fc GunnlllQ^ 




1^22 1 JdrgjLgji_S:l,-U^QQ;'^ lYang et al.l 1^04 iRIake et al.l 

2004 ; lTVaiietaLn200^ have used [O II] as the marker, 
either because it is the only emission line diagnostic 
available at high redshifts or for the sake of facilitating 
comparisons between low and high redshifts and among 
different authors. Recently, several groups llGoto et al.l 
1200.4 [Quintero et al.ll2004 IBalogh et al.ll200.^ have em¬ 


ployed Ha emission as the marker for star formation, 
rather than [O II], for samples of galaxies from the Sloan 
Digital Sky Survey (SDSS). Usually, a low equivalent 
width (EW) threshold on ]0 II], e.g., 2.5A, is used as 
the criterion for non-detection of such emission lines, ac¬ 
cording to the sensitivity in each survey. This reflects an 
underlying assumption that any emission line would be 
coming from starforming regions. If this assumption fails 
to hold—for instance, if the contribution of [O II] from 
AGN is significant—the post-starburst sample dehned by 
such a method would be incomplete because AGN would 
be mistakenly discarded as star-forming galaxies. Thus, 
using ]0 II] in post-starburst galaxy studies requires a 
better understanding of the many origins of JO II] in all 
types of galaxies. 

In many studies of AGNs, the opposite question is 
asked: what fraction of ]0 II] emission in AGN spec¬ 
tra comes from star formation? A variety of studies 
of high-ionization AGN (QSOs and Seyferts) have at¬ 
tempted to separate the contributions from AGN and 
star formation to line emission. It is of special interest 
to the understanding of t he coevolution of th e nuclear 
black holes and the bulges. iGroom et al l ll2002ll found by 
measuring the Baldwin effects llBaidwmltl977ll that most 
[O II] emission in QS Os might be from the h ost galaxies 
instead of the AGN. iRichards et al.l ll200,4l found that 
dust-reddened QSOs have stronger JO II] emission than 
normal QSOs. One of the many possible explanations 
is that the host galaxies of those QSOs with higher ex¬ 
tinctions have higher star formation rates. In contrast, 
studies by |h 3 il2nnl show that in quasar spectra there 
is very little [O II] emission beyond that expected from 
the AGN itself, indicating a supressed star formation ef¬ 
ficiency in quasar host galaxies. 

Unlike the controversy in Type I QSOs, the picture 
is a little clear er for Seyfert 2s and Type II QSOs. 
iGii et all i200(lll found that [O IIIJ/H/3 and [O III]/[0 II] 
ratios are lower in those Seyfert 2 galaxies with signif¬ 
icant star formation, which indicates that star forma¬ 
tion could contribute substantially or even dominate the 
H/3 and [O II] emission. Gonsistently, iKim et a T] il2nn3i 
found that Type II QSOs exhibit significantly enhanced 
[O II]/JO III] ratios relative to Type I QSOs. Combined 
with their high ]0 III] luminosities, the high [O II]/JO III] 
ratio is best explained by a high level of star formation. 
However, the problem is far from settled. Although the 
[O II] emission in QSO spectra could be a combination of 
two or many origins, we still have no idea of their propor¬ 
tions. In addition, QSO hosts are rare among galaxies. 
Narrow-line, low-luminosity AGN and LINERs (which 
may or may not be AGN) are much more abundant. The 


dominant origin of ]0 II] emission among these galaxies 
is still unknown. 

A variety of approaches can help to distinguish emis¬ 
sion lines having different origins. One commonly used 
method, w hich is relatively re l iable, is emission-line 
diagnostics llBaldwin et alJIlQSH IVeilleux & OsterbrockI 
Il987t iRola et al.lll997l: iKauffmann et al.ll200.4 etc.b as 
also used by many authors mentioned above in QSO 
studies. This method is effective in distinguishing the 
two major origins of emission: star formation versus 
AGN. We will focus on using this method to investigate 
the origin of [O II] emission in red galaxies in this paper. 

The large SDSS redshift survey provides a perfect sam¬ 
ple for such a study. Its wide wavelength coverage covers 
many emission lines from [O II] to Ha, from which many 
line ratio diagnostics can be created. At the same time, 
moderate resolution still allows reasonable line profile fit¬ 
ting to be conducted, giving good control of errors in the 
line strength measurement. Finally, its unprecedented 
huge sample size produces reliable st atistics. Previous 
studies of line ratio diagnostics fe.g.. IRola et al.iri997ll 
used much smaller datasets with less uniform sampling. 

In this paper, we present the discovery of a bimodality 
in [O II]/Ha ratio among galaxies. One mode is largely 
associated with star-forming galaxies. Galaxies in the 
other mode generally have line ratios similar to LINERs, 
with ]0 II] emission produced by a mechanism not associ¬ 
ated with ordinary star formation. Narrow-line Seyferts 
and Transition Objects mostly fall in between the two 
dominant populations, consistent with the picture that 
both star formation and AGN (or some other source) 
might contribute substantially to the emission in these 
objects. The [O II]/Ha bimodality we find also has im¬ 
portant implications for post-starburst studies. 

The paper is structured as follows. In we describe 
the data used. (The details of our measurement methods 
are descibed in the Appendices.) In we compare [O II] 
emission-line strength with Ha for all types of galaxies. 
We demonstrate that red and blue galaxies show different 
[O II]-Ha correlations, which appear to reflect different 
emission origins. Glassification based on this ] 0 II] /Ha 
bimodality are made and investigated later in iI4.ll with 
the standard line ratio diagnostics. Possible origins of the 
emission in red galaxies are discussed. With the conclu¬ 
sions drawn from that, we show in Jj^lhow the selection of 
post-starbursts using [O II] leads to an incomplete sam¬ 
ple and explore the nature of emission in post-starbursts. 
We summarize in The cosmology used is a flat A cold 
dark matter (AGDM) cosmology with density parameter 
nm = 0.3. 

2. OVERVIEW OF DATA AND LINE MEASUREMENTS 

The SDSS llYork et al.ll200nl: iStoughton et aljr2002ll is 
an ongoing imaging and spectroscopic survey that will 
eventually cover ^ tt steradians of the celestial sphere. 
It utilizes a dedicated 2.5-m telescope at Apache Point 
Observatory. The imaging is done with five broadband 

g, r, i, and z; iFukugita et alJ 
. Spectra are obtained with 
two fiber-fed spectrographs, covering the wavelength 
range of 3800-9200A with a resolution of R ~ 2000. 
Fibers have a fixed aperture of 3”. 

The spectroscopic data used here have been re¬ 
duced through the Princeton spectroscopic reduction 


filters in drift scan mode (u , 
' iton et aTIl 
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pipeline (Schlegel et al. in prep), which produces 
the flux- and wavelength-calibrated spectra^. The 
redshift catalog of galaxies use d is from th e NYU 
Value Added Galaxy Gatalog llBlanton et al.l l2005bL 
http: //wassup.phy sics.nYU.edu/vagc/). K-corrections 
were derived using iRlanton et all ll2003aD ’s kcorrect code 
v3_2. All the magnitudes used in this paper are K- 
corrected to a redshift of z = 0.1 by shifting the 
SPSS bandpasses to b luer wavelengths by a factor of 1.1 
llBlanton et al.ll2003bf) . Thus all magnitudes and color 
notations carry a superscript or subscript of 0.1 . 

The results in this paper are based on the spectra of 
~ 40 0,000 galaxies contained in the S PSS Pata Release 
Four I Adelman-McGarthv et alJ2000l PR4). We employ 
a catalog of objects targeted by the main galaxy survey 
that have been spectroscopically confirmed as galaxies; 
QSOs are not included in this sample. T arget selection 
for th e main galaxy sample is described in iStrauss et al.l 
(I 200 I . The magnitude limit is r = 17.77 in Petrosian 
magnitudes after c orrection for the foreg round galactic 
extinction following iSchlegel et alJ ill 9981) . 

For many plots (Figs. 12 0 0 EJ El and El in this 
paper, we used a volume-limited sample of about 55,000 
galaxies with 0.07 < z < 0.10 and Mo. 1,. —5 log h brighter 
than -19.5, and to have well-measured spectra in the 
vicinity of [O II] and Ha (using criteria described in the 
Appendix 0j). However, in color-magnitude diagrams 
and line ratio diagnostic plots (Figs. 00 [30 0^2 and 
d, a magnitude limited sample with 0.05 < z < 0.10 
and r brighter than 17.77 is used to cover a wider range 
in galaxy luminosity and metallicity. The redshift ranges 
of these samples are chosen for several reasons. First, se¬ 
lecting very nearby galaxies is unfavorable because Sloan 
spectra are taken with a 3” fiber which only covers a 
small fraction of each galaxy at low redshifts. Secondly, 
objects at redshifts greater than z = 0.15 are disfavored 
because blue galaxies drop out the r-band selected sam¬ 
ple quickly with increasing z. Additionally, the bad col¬ 
umn on the CCP around 7300A heavily hinders the accu¬ 
rate measurements of Ha at z ^ 0.114 and other lines at 
higher z. The median S/N of the spectra in the volume- 
limited sample is 18.5 per pixel. 

For each galaxy, we subtract the stellar continuum and 
measure the emission line fluxes and equivalent widths 
for [O II] A3727, H/3, [O III] A5007, [O I] A6300, Ha, and 
[N II] A6583. Zero point errors in the equivalent widths 
are removed and the corrections are propagated into 
fluxes. Reliable emission line detections are defined as 
3 or more a detection in EW. For the technical details 
of our stellar continuum subtraction procedure, emission 
line measurements, and zero point determinations, see 
Appendix 0] and 0| 

3. [O II] A3727 VS. Ha 
3.1. Color Bimodality 

Before diving into line comparisons, we recall that 
the color distribution of galaxies i s highly bimodal. as 
have been s hown r epeatedly (e.g .. iStrateva et al.l 120011: 
Hogg e^^ _|2002_|Blanton_e£ay_j200^ iBaMr^ 

H Bell et al.ll2004 IWeiner et alJl2005t Iwillmer et al. 

Giallongo et al.ll2005D. Fig. 0 shows the color- 

magnitude diagram for all SDSS galaxies within 0.05 < 



Fig. 1.— Rest-frame color-magnitude diagram for all SDSS 
galaxies with redshifts between 0.05 and 0.10. The two solid lines 
indicate the cut we adopt for separating red and blue galaxies. Red 
galaxies are defined as the points above the upper line, while blue 
galaxies are below the lower line. We leave a gap in between to 
eliminate galaxies with potentially ambiguous status. The vertical 
dashed line indicates the absolute magnitude limit for the volume- 
limited sample. 


z < 0.10. We use this sample rather than the volume- 
limited sample simply to display a wider range in galaxy 
luminosity. 

In color-magnitude space, galaxies fall mainly into two 
categories. One category is the tight sequence between 
color and magnitude at the red end (larger value of 
— r)), commonly called the “red sequence”; the 
other is the swath of points at bluer colors, which is 
sometimes referred to as the “blue cloud”. Empirically, 
blue galaxies mostly have disk-dominated morphologies 
and are actively forming stars; while red galaxies mostly 
have bulge-dominated morphologies and are relatively 
quiescient in terms of star formation. 

Although the bimodality is obvious, it is impossible to 
entirely separate the two populations by color. One pos¬ 
sibility is that color provides a poor measurement of some 
more fundamental property that cleanly distinguishes 
galaxies into two distinct populations. Alternatively, the 
absence of a clean separation in color may simply reflect 
the smooth evolution of one population into the other. 
It is also likely that both possibilities are involved. For 
now, we use conservative color cuts to define the samples 
of red and blue galaxies, as shown with the two solid 
lines in Fig. 0 and desribed by the following inequali¬ 
ties. We leave a gap in color between the two samples to 
reduce the importance of color errors and galaxies with 
ambiguous status. We define these color cuts by: 

Red: °-^(g - r) >-0.025(Mo.i^ - 5 log h)-H 0.42 (1) 
Blue : - r) < -0.025(Mo.i,. - 5 log h) + 0.35. (2) 


^ http://spectro.princeton.edu/ 
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3.2. [O II]/Ha Bimodality 

In this paper, we will use both equivalent width (EW) 
and flux/luminosity measurements of emission lines. For 
line intensity ratios and SFR measurements, the use of 
flux is necessary, but fluxes are much more sensitive to 
dust extinction than equivalent widths^. Furthermore, 
EW will not be affected by imperfect flux calibration. On 
the other hand, to interpret the EW ratio between two 
lines, the continuum shape (i.e., the color of the galaxy) 
needs to be taken into account. In addition, line luminos¬ 
ity and EW describe different sorts of physical quantities: 
line luminosity describes the total amount of emission, 
while EW reflects the strength of emission compared to 
the luminosity of the galaxy. We will primarily employ 
EW in this paper because it is not subject to variations 
in galaxy size (in contrast, line luminosity tends to scale 
with the total luminosity of a galaxy). It is worth noting 
that the correlations in EW are often tighter than those 
in luminosity, for all these reasons. 

Figure 121 compares [Oil] EW with Ha EW for galax¬ 
ies in the volume-limited sample (0.07 < z < 0.10, 
Mo.ir — 51og/i < —19.5) which have both [O II] and Ha 
positively detected. Clearly, there are two [O II]-Ha se¬ 
quences. One has a shallow slope and extends to very 
high Ha EW. The other sequence has a rather steep 
slope and very high concentration of points at low Ha 
EW, with [O II] EW spanning a wide range. 

The bimodality in the [O II]-Ha diagram becomes even 
more obvious if plotted on a logarithmic scale, as shown 
in the upper right panel of Fig.|21 Here the shading re¬ 
flects the density of points in units of number per square 
dex. The Ha EW distribution can be well fit by two log¬ 
normal distributions, one centered around lA and the 
other centered around IdA. In linear space, the high-EW 
population spans a wide range of EWs, but it becomes 
much more concentrated in log space. 

The bottom two panels in Fig.[2|show [O II] EW vs. Ha 
EW for the red galaxies and blue galaxies separately. The 
clear difference between the two populations shows im¬ 
mediately that most red galaxies reside on the steep-slope 
sequence, while most blue galaxies are on the shallow- 
slope sequence. Clearly, the bimodality in [O II]/Ha ra¬ 
tio echoes the color bimodality. Since most blue galax¬ 
ies are actively forming stars, this is telling us that the 
shallow-slope sequence must have an [O II]/Ha ratio con¬ 
sistent with that produced in HII regions photoionized by 
O and B stars in star-forming galaxies. What mechanism 
gives rise to the different ratio seen in red galaxies? 

Before addressing this question, we detour briefly to 
consider another issue. Since we are thus far using 
equivalent width instead of line luminosity, one might 
imagine that the bimodality has little to do with emis¬ 
sion strength; rather, differences in eqivalent-width ratio 
might reflect only the different continuum ratios between 
stellar continua near [O II] and Ha in red galaxies ver¬ 
sus that of blue galaxies. Since red galaxies have higher 
continuum ratios of E’c(Ha)/Fc([0 II]), they will have a 
higher value of EW([0 II])/EW(Ha) for the same emis¬ 
sion line luminosity ratio. This possibility can be checked 
by comparing L([0 II]) to L(Ha) directly. 

^ This is not necessarily true in cases when the emission line 
region and the stars dominating the continuum have significantly 
different dust obscuration. 


As shown bv iKennicuttI iHQflflfl . however, such a plot 
of luminosity versus luminosity can be quite misleading 
because of the size effect: for a given class of galaxies, 
we expect both L([0 II]) and L(Ha) to each have a sig¬ 
nificant correlation with the stellar mass of the galaxies, 
especially at the low-mass end, which would veil the true 
relationship between [O II] and Ha. This hidden variable 
can be easily removed by using specific line luminosity; 
here we divide the line luminosity by the Sloan band 
luminosity, which is a coarse proxy for the stellar mass. 
This effectively creates an EW, but one that does not 
include the continuum color variation between red and 
blue galaxies. As shown in Fig. EKa), the same sepa¬ 
ration between the two sequences persists in a plot of 
L([0 II])/L(°-^r) vs. L(Ha)/L(°-^r), but with a smaller 
opening angle between them. The wider opening angle 
in the equivalent width plot indeed reflects the different 
continuum ratios between blue and red galaxies, but this 
is not the whole origin of the bimodality. In panel (b) 
of Fig. 121 we plot L([0 II])/L(°-^r) vs. L(Ha)/L(°-^r) in 
logarithmic scale, analogous to the upper right panel in 
Fig.m The bimodality in [O II]/Ha ratio again stands 
out clearly. We omit the plots for subsamples separated 
in color, due to their close resemblance to the bottom 
panels of Fig. [21 

Since the continuum ratios are not responsible for the 
[O II]/Ha ratio bimodality, the question of the source of 
line emission in red galaxies persists. Can the difference 
between line ratios for blue and red galaxies in [O II]/Ha 
ratio arise from reddening or metallicity dependences for 
star-forming regions? We investigated this possibility us- 
ing the SFR(Hc^ and SFR([0 II]) calibrat ions derived by 
iKennicutil lll998r and iKewlev et al.l ll2004ll , which explic¬ 
itly take into account extinction and metallicity correc¬ 
tions: 


SFR(Ha) = 7.9 x 10-'^^L(Ha/eTgs-/MQyr-^ (3) 


and 


SFR([0 II]) = 


7.9x10 4^L([0 II]/ergs / , 

16.73- 1.75[log(0/H)-t 12]^®^ 


By identifying SFR(Ha) with SFR([0 II]), we get 


L([0 II])/L(Ha) = 16.73 - 1.75[log(0/H) -f 12]. (5) 


First, we can check that the [O II]/Ha ratio seen among 
blue galaxies is consistent with that produced in star¬ 
forming HII regions. The galaxies to the right and below 
the demarcation in Fig. |21have a median [O II]/Ha flux 
ratio of 0.317, as indicated by the lower dashed line in 
Fig. 121). Their median Ha/H/3 ratio (the Balmer decre¬ 
ment) is 4.75. This corresponds to a median extinction 
Av of 1.60 (assuming Rv = Av/E{B — V) = 3.1 and 
an intrinsic Ha/H/3 ratio of 2.85 for case B recombina- 
tion at T — lO^K and Ue ^ 10^ — 10^cm“^ llOsterbrockI 
(ilia). Applying the resulting extinction correction, we 
obtain a median intrinsic [O II]/Ha ratio of 0.918. Ac¬ 
cording to Eq.|21 this [O II]/Ha ratio corresponds to a gas 
phase metallicity (log(0/H) -I- 12) of 9.0 3, which is typi¬ 
cal am ong star-forming galaxies in SDSS llTremonti et alJ 

l2?iril . 

Next, we investigate the [O II]/Ha ratios among the 
galaxies to the left and above the demarcation in Fig. |21 
These galaxies have a median [O II]/Ha flux ratio of 
1.38, as indicated by the upper dashed line in Fig. I21)- 
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Fig. 2.— (a) The distribution of all SDSS galaxies in the volume-limited sample with 0.07 < 2 < 0.10 in Hck and [O II] equivalent width 
for those objects where both lines are detected. There are two [O Ilj-Ha sequences with different slopes. The dashed line indicates an 
empirical cut (Eq|^ separating the two sequences and is reproduced in following panels. We refer to the population to the left of the line as 
“High-[0 II]/Hq: galaxies”, and refer to the population to the right as “Low-[0 Ilj/Ha galaxies”, (b) Same as panel (a) but in logarithmic 
space. The shade indicates the number of points per sq. dex. The bimodality in [O II]/Ha ratio is even more obvious in this panel. The 
dashed curve corresponds to the dashed line in panel (a), (c) Same as (b) but only for red galaxies, most of which are on the steep-slope 
sequence. Shade scale is the same as in panel (b). (d) Same as (b) but only for blue galaxies, most of which are on the shallow-slope 
sequence. These generally have an [O II]/Ha ratio consistent with an origin from star-forming HII regions (see id.21 . Dashed lines in all 
four panels indicate the same line, which we use throughout this paper to separate the two populations. 


The median Ha/H/3 ratio for these objects is 4.46, 
which corresponds to a median extinction in Ay of 
1.40 (following the same assumptions as for star-forming 
galaxies). Applying the resulting extinction correction, 
we obtain a median intrinsic [O II]/Ha ratio of 3.54. 


However, iKewlev et al.l ll2004^ shows that the maxi¬ 
mum [O II]/Ha ratio which can be reached by star¬ 
forming regions is around 2.1, far below the ratio we 
find among red galaxies. If we were to take into account 
the metallicities inferred from the [N II]/[O II] ratio or 












6 


Yan et al. 


n 

O 

I 

X 


u 



10 
8 
6 
4 
2 
0 

-2 

-5 0 5 10 15 20 

L(Ha)/L(°-‘r) xlO^ 



100.0 


10.0 


1.0 



Fig. 3.— (a) [Oil] line luminosity vs. Ha line luminosity for all galaxies in the volume-limited sample with 0.07 < z < 0.1 detected in 
both lines. Line luminosities are divided by the SDSS band luminosity to take out the size effect (see E3- The same bimodality 
exists as seen in EW plot, (b) Same as panel (a) but in log-log space. The bimodality clearly stands out. The shade indicates the number 
of galaxies per sq. dex. The two dashed lines show [O II]/Ha ratios of 1.38 and 0.317, corresponding to the median ratio for each of the 
two modes, respectively. The effect of Ay = 2 of extinction (for Ry = 3.1) is indicated by the arrow, as we do not correct for reddening. 


([O II] A3727+[0 III] AA4959,5007)/H/3 {R 23 ) among the 
red galaxies, the discrepancy would become even larger. 
As we will show in using full line ratio diagnostics, 
red galaxies in fact have line ratios characteristic of LIN- 
ERs/AGNs, not star formation. Hints of this possible 
AGN origin of high [O II]/Ha rat ios can also be seen in 
earlier work bv iRola et alJ l|1997jl . who suggested a new 
line diagnostic method using EW(H/3)and EW([0 II]). 

In Fig. 121 a demarcation is shown separating the two 
[O II]-Ha sequences; it is defined by 

EW([0 II]) = 5EW(Ha) - 7. (6) 

There is some overflow of points across the demarcation 
in the bottom panels of Fig.|21 Roughly 3% of blue galax¬ 
ies in the volume-limited sample sit to the left of the de¬ 
marcation. Nearly half of these are post-starburst galax¬ 
ies, which have had their star formation quenched to a 
fairly low level but still look blue because they contain 
relatively young stars. For the red galaxies, a much larger 
number of galaxies(^ 30%) lie to the right of the line. 
As will be shown in sm based on their other line ra¬ 
tios, most of these galaxies are Transition Objects (TO), 
while a small fraction are Seyferts and dusty starforming 
galaxies. They will also be investigated further below. 

3.3. Classification by [O II]-Ha 

So far, we have been studying the galaxy distribution 
in [O II]-Ha EW plot for subsamples separated in their 
rest-frame colors. While in selecting our blue and red 
subsamples we exclude a gap in — r) color to reduce 
contamination, the overflow of points across the [O II]- 
Ha demarcation is still significant for red galaxies. We 
can also investigate the bimodality from the other direc¬ 
tion by dividing the sample according to position on the 
[O II]-Ha EW plot and examine the distribution of each 
subsample in color-magnitude space. 


We show the demarcation defined above in Eq. El as 
the dashed line in Fig. |21 From now on, we refer to 
the population to the left of the line as “High-[0 II]/Ha 
galaxies”, and refer to the population to the right as 
“Low-[O II]/Ha galaxies”. Note that there is a large 
population of galaxies which do not have both [O II] and 
Ha detected and are not plotted here. Their EW mea¬ 
surements, if intrinsically zero, should have a symmetric 
Gaussian distribution around zero. These galaxies can be 
again classified into three groups based on the emission 
lines detectable in them: 

1. Galaxies with [O II] undetected and Ha detected. 
We label these “Ha-only galaxies”. Note that other 
emission lines may be present in these objects’ 
spectra besides [O II]. Similarly, the “[O II]-only” 
galaxies defined below may have detectable lines 
besides [O II], but have no significant Ha emission. 

2. Galaxies with [O II] detected and Ha undetected. 
These galaxies likely belong to the same population 
as the “High-[0 II]/Ha” galaxies and fall to the left 
of the demarcation line. We label these “[O Ilj-only 
galaxies”. 

3. Galaxies with neither [O II] nor Ha detected. We 
label these galaxies as “Quiescent Galaxies” in the 
remainder of this paper. Their EW measurements 
form a symmetric distribution around (0,0) in the 
[O II]-Ha diagram and fall to the left of our de¬ 
marcation, as shown in the inset plot of Fig. H As 
illustrated in Fig. 01 quiescent galaxies have a color- 
magnitude distribution very similar to the “High- 
[O II]/Ha” galaxies. 

Fig.EI shows the color-magnitude diagrams for these 
subsamples: quiescent galaxies (panel a), the high- 
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[O II]/Ha galaxies (b), [O II]-only galaxies (c), the low- 
[O Il]/Ha galaxies (e), and Ha-only galaxies (f). Our 
color cuts are drawn here as the two straight lines for 
references. Both the quiescent galaxies and the high- 
[O II]/Ha galaxies display a highly symmetric red se¬ 
quence. The quiescent galaxies have a wider luminos¬ 
ity range and the high-[0 II]/Ha galaxies have a slight 
asymmetry in color, but they are otherwise nearly in¬ 
distinguishable in color-magnitude space. As shown in 
Fig-El they also have similar Sersic index ^ distributions. 
As expected, the [O II]-only galaxies have similar color- 
magnitude distributions and Sersic index distributions 
as the high-[0 II]/Ha galaxies. The combination of the 
three subsamples shown in Fig.^Ji-c - i.e., those galaxies 
to the left of the demarcation given by Eq.Elon [O II]-Ha 
diagram - together comprise a tight red sequence in the 
color-magnitude-concentration space, as shown in panel 
(d) of Fig. 0] 

The low [O II]/Ha subsample, as expected, covers the 
blue cloud in the color-magnitude diagram (Fig. 01i). 
However, it also includes some galaxies with red- 
sequence-like colors, especially at the faint end. These 
red galaxies within the low [O II]/Ha sample, which in 
fact are the overflow points in the lower left panel of 
Fig.m also deserve special investigation. We will re¬ 
fer to them as “Red-Low-[O II]/Ha galaxies”. As shown 
in Fig. 0 these galaxies generally have lower Sersic in¬ 
dices than the overall red galaxy population. Their dis¬ 
tribution does not resemble either of the possible parent 
samples; i.e., the red galaxy sample or the low-[0 II]/Ha 
sample. Thus, although these galaxies have red colors, 
they are likely to be a different population from the red 
sequence. The Ha-only galaxies (panel f) also extend 
to red colors. Without detection of [O II], their catego¬ 
rization is difficult; we will discuss them briefly later on. 
To recapitulate, we list in Tabled our sample definitions 
and the fraction of each category in the volume-limited 
sample. 

The [O II] EW distribution of all the galaxies to the 
left of the demarcation can be modeled well by the com¬ 
bination of a Gaussian distribution centered near zero 
and a log-normal distribution, as shown in Fig. The 
two components correspond closely to the EW distribu¬ 
tions of the quiescent population alone and the com¬ 
bined population of the high-[0 II]/Ha and the [O II]- 
only galaxies, respectively; the fractional area covered by 
the log-normal distribution (35.4%) is nearly identical to 
the fraction of galaxies with detectable [O II] emission 
(35.5%). 

In the following section, we employ more line ratio di¬ 
agnostics to investigate the source of the emission lines in 
red galaxies, specifically, the High-[0 II]/Ha population 
and the Red-Low-[O II]/Ha population. 


4. ORIGIN OF [Oil] EMISSION IN RED GALAXIES 
4.1. Classifications by line ratio diagnostics 


® Sersic index is a seeing-corrected structural parameter describ¬ 
ing the radial light profile of a galaxy. A Sersic index of 4 cor¬ 
responds to a de Vaucouleurs profile, and an index of 1 corre- 
ponds_to_^^n_ex£onenti^profile. The indices used here are derived 
by |Bla£ton_eL_^ J2005^. One caveat is that these Sersic indices 
are underestimates at high Sersic index, e.g., n=3.5 for a de Vau¬ 
couleurs profile. 


Emission lines in galaxy spectra can be produced by 
a variety of mechanisms; Among them, photoionization 
by AGN and by hot stars in star-forming HII regions 
are tw o maj or mechanisms. Baldwin, Phillips and Ter- 
levich ill 98lL hereafter HPT) pioneered the use of inten¬ 
sity ratios of two pairs of emission lines to distinguish 
starforming galaxie s from AGNs. The me t hod w as fur¬ 
ther developed by iVeilleux fc OsterbrockI lll987ll . The 
so-called HPT diagram of line ratios, such as the one 
shown in the panel (a) of Fig. 0 has become a standard 
tool for determining the source of line emission. 

All the galaxies plotted in Fig.dhave at least 2cr detec¬ 
tions of all emission lines considered (H/3, [O III] A5007, 
Ha, [N II] A6583). Two sequences of emission-line galax¬ 
ies dominate this diagram; one corresponds to the star¬ 
forming galaxies, and the other to AGN. The star¬ 
forming sequence is the curved concentration of points 
that extends from center bottom to upper left; its shape 
reflects the correlation between ioni zation parameter an d 
metallicity in star-forming galaxies ijKewlev et al.ll200Hi . 
The concentration of points on t he right part of the 
plot is identified with AGNs fe.g.. iBaldwin et alJlT 981l: 
IVeilleux fc Osterbrocklll987HKauffmann et al.ll2003H. In 

studies employing optical emission lines from large aper¬ 
ture spectroscopy, one usually cannot confirm that the 
emission lines are associated with a central supermassive 
black hole. From now on, we define a “true AGN” as a 
system in which the radiation is ultimately powered by 
accretion onto a supermassive black hole. In the remain¬ 
der of this paper, we will refer to galaxies with line ratios 
characteristic of true AGN as “AGN-like galaxies”. 

The theoretical modeling of starburst galaxies can pro¬ 
duce a boundary for the locus of star-forming galax¬ 
ies on the BPT diagram, which is called the extreme 
star burst classification l ine. The classification line given 
by iKewlev et alJ ll2001ll is indicated by the solid curve 
in Fig. d As seen in this figure, that would serve as 
a very conservative cut for selecting AGN-like galaxies. 
An alternative empirical dem arcation has been suggested 
by iKauffma.uu et all jl^2nn.3li , which is indicated by the 
dashed curve in Fig. |3 

Traditionally, AGNs are further divided into sev¬ 
eral sub-categories according to their positions in 
the BPT diagram: Seyferts, Low-Ionization Nuclear 
Emission-line Regions (LINERs) and Transition Ob¬ 
jects (TOs). Traditional definitions for Seyferts and 
LINE Rs are indicated by the horizontal and vertical 

lines l)Shiider V OsterbrocUl 98l( IVeilleux fc OsterbrockI 

IT^ in Fig, d Ever since the discovery of LINERs 

llHeckma,nl1 98(11 . it has remained uncertain whether they 
are true AGN. Recently, it has become clear that at 
least some LINERs are indeed accretion-powered sys¬ 
tems llFi]inneuknll2nn3l: lHnll2nn4 and references therein). 
However, LINER-type line ratios are also found com- 
ing from extended regions in some early-type galaxies 


_ re gi( _ _ 

lllleckman et al.ll989tiGoudfrooii et al.ll993:lGoudfrooil 

EST in support of other ionization mechanisms, such 

as shockwaves, cooling flows, hot old stars, etc. Tran¬ 
sition Objects are galaxies that satisfy most of the 
LINER criteria but have weaker [O I] A6300 /Ha ratio 
ifFilippenko fc Terlevichl 119^ iHo et al.lll99^ . The di¬ 
viding line between LINERs and TOs in [O I]/Ha is in¬ 
dicated by the vertical dashed line in Fig.0 which shows 
another useful BPT diagram with [O III] A5007/H/3 plot- 
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Fig. 4.— Color-magnitude diagrams for the quiescent galaxies (a), the high-[0 II]/Ha galaxies (b), the [O II]-only galaxies (c), the 
galaxies to the left of the [O II]-Hq demarcation (d), the low-[0 IIJ/Ha galaxies (e), and the Hoi-only galaxies (f) in the magnitude limited 
sample (0.05 < 2 : < 0.10). Quiescent galaxies, the high-[0 galaxies, and the [O II]-only galaxies all follow a tight red sequence. 

Low-[0 II]/Ha galaxies cover not only the blue cloud but also part of the red sequence region. Contours are plotted at number densities 
of 1.6 X lO^^n (halved in the panel b, c, and f) galaxies per mag^, where n=l,2,..., starting from outermost contour. The inset in panel (a) 
shows the [O II]-Ha EW distribution for the quiescent galaxies. They form a symmetric distribution around zero as expected. 


Name 

Nonzero? 

Criteria 

Percentage 

fo II] 

Ho 

High-[0 II]/Ha 

Y 

Y 

EW([0 II]) > 5EW(Ha) - 7 

13.3% 

Low-[0 II]/Ha 

Red-Low -[0 II]/Ha 

Y 

Y 

EW A H]) < 5EW(Ha) - 7 

(as low-[0 II]/Ha, plus 
-r)> -0.025(Mo.i^ - 5logh) -|-0.42) 

38.0% 

(4.4%) 

[0 II]-only galaxies 

Y 

N 


5.2% 

Ha-only galaxies 

N 

Y 


15.7% 

Quiescent galaxies 

N 

N 


27.9% 


TABLE 1 

Subsample definitions and their fractions in the volume-limited sample; see El FOR DETAILS. 


ted against [O I] A6300/Ha. The nature of Transition 
Objects still remains unclear. 

Now we examine the position of the two classes of red 
galaxies in the BPT diagrams. We look at the High- 


[O II]/Ha sample and the Red-Low-[O II]/Ha sample 
separately because they have different [O II]-Ha corre¬ 
lations. In the two BPT diagrams (Fig- EI&IHIlj the two 
samples are displayed by themselves in panels (b) and 
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Fig. 5. — Sersic index distributions of galaxies for four sub¬ 
samples as defined in Table ^ high-[0 II]/Ho galaxies (thick- 
solid, height doubled for easy comparison), low-[0 II]/Ha galax¬ 
ies (thick dashed), quiescent galaxies (thin-solid), and red-low- 
[O II]/Ha galaxies (thin dashed). The high-[0 II]/Ha galaxies 
have a similar distribution to the quiescent galaxies, and very differ¬ 
ent from that of the low-[0 II]/Ha galaxies. The red-low-[0 II]/Ha 
galaxies have a distribution unlike either possible parent sample. 


(c). In both BPT diagrams, the High-[0 II]/Ha sample 
(panel b) has a nearly symmetric distribution centered 
on the region classified as LINERs. The scatter is due to 
intrinsic variations as well as measurement errors. The 
median errors are indicated on the plot. Galaxies in the 
Red-Low-[O II]/Ha sample mostly fall into the transition 
region in both BPT diagrams (panel (c) in each). Fig.|Hl 
shows clearly that most of the se galaxie s are in f act Tran¬ 
sition Objects as defined bv iHo et~ ] lll997ar . while a 
smaller fraction are Seyferts and star-forming galaxies. 

Gomparing [O II] A3727 with [O I] A6300, one can see 
the similarity in their relations to Ha: both [O II]/Ha 
and [O I]/Ha ratios provide a separation between LIN¬ 
ERs and other categories. In fact, the bimodality in 
[O I]/Ha is employed in one of the BPT diagrams to 
distinguish AGNs from star-forming HII regions. Un¬ 
fortunately, [O I] A6300 is often too weak to detect, 
which hinders the use of this particular BPT diagram 
for the full sample. As we have shown above, the 
strong bimodality in EW([0 II])/EW(Ha) (or similarly, 
EW([0 II])/EW(H/3) ) might make it a good alternative, 
despite the wide separation of the two lines in wave¬ 
length. We propose the use of the demarcation in the 
EW([0 II])-EW(Ha) plane shown in Fig. |21 to separate 
LINER-like galaxies from other line-emitting objects. 

Not all of the high-[0 II]/Ha galaxies are plotted on 
the two BPT diagrams(Fig.[7|&|3). Because [O I] A6300 
and H/3 are much weaker than strong emission lines 
like [O II] A3727 and Ha, requiring them to be detected 
significantly reduces the sample size, especially among 



0 pdf I _^_1_L.^_1_^ I _1_ I .. . 

-5 0 5 10 15 20 

[on] EW (A) 


Fig. 6.— The histogram shows the [O II] A3727 EW distribution 
of all galaxies to the left of the demarcation (Eq. in [O II]- 
Ha diagram, including galaxies consistent with zero EW. It can 
be modeled well by the combination of a Gaussian distribution 
centered near zero (the dashed curve) and a log-normal distribution 
(the dot-dash curve). The sum of the two components is shown 
as the solid curve, which goes through the histogram. The EW 
distributions of all other emission lines in this subset of galaxies 
can be modeled the same way, with similar results. As in the rest 
of this paper,plotted EWs are the values after zero point correction 
(see Appendix I bI for details). 


the high-[0 II]/Ha population. Do such requirements 
bias our sample selection? We test this with a BPT 
diagram (Fig. |3 which only involves the strong lines: 
[OlII]/[OlI] vs. [N II]/Ha. Panel (a) of Fig. 0 shows 
the distribution for all galaxies with at least 2-cr detec¬ 
tions on all of the four strong lines. Galaxies identified as 
star-forming from other BPT diagrams are again found 
from bottom to middle left. The horizontal lin e indicates 
the [O III]/[0 II] division originally defined bv lHeckmanl 
iTT^ to separate LINERs from Seyferts. The vertical 
line is a commonly adopted limit for separatin g LINERs 
and S eyferts from star-forming galaxies (e.e.. IHo~ iil] 
ll997aH . In panel (b) of Fig. El only high-[0 II]/Ha galax¬ 
ies are plotted. It is worth noting that the high-[0 II]/Ha 
galaxies again fall into the LINER region. In panel (c), 
we only plot high-[0 II]/Ha galaxies which also have at 
least 2-(T detections on Jp I] and H/3, i.e., those which 
are also plotted in Fig. [7|&|8| These galaxies have very 
similar distributions as those without detectable [O I] or 
H/3, only biased slightly against high [N II]/Ha ratio and 
high [O III]/[0 II] ratio. The former is likely due to the 
difficulty of measuring H/3 for galaxies with very low Ha 
and thus high [N II]/Ha ratio. The latter is probably due 
to the difficulty of measuring [O I] for galaxies with rela¬ 
tively high ionization. Thus, requiring [O I] and H/3 de¬ 
tections only biases the sample slightly, and in expected 
ways. We thus believe that nearly all high-[0 II]/Ha 
galaxies have LINER-like line ratios, regardless of their 
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Fig. 7.— (a) Line ratio diagnostic diagram (BPT diagram) showing the emission-line flux ratio [O versus the ratio [N II]/Hq 

for the magnitude-limited sample (0.05 < 2 : < 0.1). This BPT digram, by construction, is very insensitive to the effects of e xtinction. All 

S ’3s plotted have at least 2g detection on all four emission lines used. The solid curve shows th e demarcation def i ned b y lKewlev et ^ 
to separate star-forming galaxies from AGN. The dashed curve shows the demarca tion use d bylKauffinamLA-aU I2QQ3) • Conventional 
divid ing; lines for Seyferts and LINERs are indicated by the horizontal and vertical lines iShuder^rOsterbroci^^^nWhleinr^ OsterbrocU 
[19^ . Contours are plotted at number densities of 24,000n galaxies per unit area, where n=l,2,..., starting from the outermost contour, 
(b) Same diagram as panel (a) but only High-[0 II]/Ha galaxies are shown. These are centered in the LINER region. The Icr measurement 
errors propagated to log-line-ratio is shown as the cross. Contours are plotted at number densities of 8, OOOn galaxies per unit area, where 
n=l,2,... . (c) Same diagram for the Red-Low-[0 II]/Ha sample. They lie to the left of the LINER region and extend down to the 
star-forming region. Other than the Seyferts and the star-forming galaxies, most are Transition Objects, as shown in Fig. 0 Contours 
indicate the same number density levels as in panel (b). 



-2.5-2.0-1.5-1.0-0.5 0.0 -( 255 - 2 . 0 - 1 . 5 - 1 . 0 - 0.5 0.0 -C255-2.0 - 1.5 - 1.0-0.5 0.0 0.5 
Log([0l] A6300/Ha) Log([0l] A6300/Ha) Log([0l] A6300/Ha) 


Fig. 8.— (a) BPT diagram showing the emission-line flux ratio [O III] A5007/H;9 versus the ratio [O I] A6300/Ha for the magnitude- 
limited sample (0.05 < 2 : < 0.1). The effect of Ay = 2 of extinction (for Ry = 3.l) is indicated by the littl e arrow. All galaxies plotted here 
have at least 2cr detections on all four emission lines. The solid curve shows the demarcation defined bv IKewlev et a y ^2o^!^ to separate 
star-forming galaxies (to the left and below) fr om AGN (to the right and ab ove) . The horizontal line is the di viding line between Seyferts 
(above) and LINER/TOs (below) as defined bv IShuder fc Usterbrocfl 1198111 and IVeilleux fc OsterbrocU 1198711 . Contours indicate number 
density levels of 10, OOOn (n=l,2,...) galaxies per unit area, (b) Th e same di a gram f or the High-[0 II]/Ha sample. The vertical line is 
the dividing line between Transition Objects and LINERs given bv lHo et all il997all . The High-[0 Il]/Ha sample lies almost entirely 
in the LINER region, with scatter similar to the measurement errors. The la measurement errors propagated to log-line-ratio is shown 
as the cross. Contours indicate number density levels of 2,400n (n = 1,2,...) galaxies per unit area, (c) The same diagram for the 
Red-Low-[0 II]/Ha sample. These objects lie to the left of the TO/LINER dividing line and thus are indeed Transition Objects, plus a 
few Seyferts and star-forming galaxies. Contours indicate the same number density levels as in panel (b). 
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[O I] or H/3 detectability. 

The high-[0 II]/Ha galaxies, as shown above, have 
relatively uniform line ratios in [O II]/Ha, [O III]/H/3, 
[N II]/Ha, [O I]/Ha, and [O in]/[0 II], w hich a ltogether 
identify them as LINER-like. As found in a the high- 
[O II]/Ha galaxies occupy the same region in the color- 
magnitude space as quiescent galaxies do and they share 
the same distribution in galactic concentration. Since 
they only differ in emission line strengths, we speculate 
that these LINER-like galaxies have a close evolution¬ 
ary relationship with the quiescent red sequence galaxies. 
They might be the immediate progenitors of the quies¬ 
cent red galaxies — a scenario in which the LINER phase 
lasts for some period of time when the galaxy first arrives 
on the red sequence (Graves et ah, in prep). Alterna¬ 
tively, it could be that LINER phases are intermittent 
through out the life of a quiescent galaxy. As pointed 
out in the combination of the LINER-like (high- 
[O II]/Ha) galaxies and the quiescent galaxies produces 
a uniform red sequence; the red sequence is effectively a 
LINER/Quiescent sequence. 

We turn now to the low-[0 II]/Ha galaxies. As was 
the case previously, not all of the low-[0 II]/Ha galaxies 
are plotted on the BPT diagrams because of S/N limi¬ 
tations. But as we found before, in the [O III]/[0 II] vs. 
[N II]/Ha diagram, the distribution of red-low-[O II]/Ha 
galaxies is also found to be independent of their [O I] 
or H/3 detectability. Their line ratios vary from star¬ 
forming, to Transition Objects, to Seyferts. To quantify 
the relative fraction of each category, we use [N II]/Ha 
ratio, which is most commonly detectable among this 
population. As listed in Tabled a little more than 1/3 of 
these have [N II]/Ha < 0.6, similar to dusty star-forming 
galaxies, while the rest have [N II]/Ha > 0.6, resembling 
TO/Seyferts. 

The [O II]-only galaxies likely belong to the high- 
[O II]/Ha category intrinsically, as they occupy the same 
color-magnitude space, have the same Sersic index distri¬ 
bution, and have lower limits on [O II]/Ha ratio higher 
than typical high-[0 II]/Ha galaxies. Therefore, we sus¬ 
pect they also belong to the LINER-like category. 

The case for the Ha-only galaxies is more complicated, 
since the ratio of the EWs at which [O II] and Ha are 
detectable is similar to the [O II]/Ha EW ratio of a typ¬ 
ical high-[0 II]/Ha galaxy. The Ha-only population will 
include both high-[0 II]/Ha and low-[0 II]/Ha galaxies. 
Among them, we can only separate star-forming galaxies 
out using the [N II]/Ha ratio, the only line ratio available 
in most of these galaxies; see Table |21 

Previous studies of emission lines in early-type galax¬ 
ies had already shown that [O II] A3727 emissio n is de¬ 
tecte d in rough l y half of all early-type gala xies llMavalll 
II9581 iGaldwelll 11984 iPhillios et alJ IT9^ . Most of 
these galaxies in which multiple e mission lines are de¬ 
tected have LINER-lik e line rati os iPhillips et al ]l1986l: 
iGoiidfrooii et alJ II994 iZeilinger et al.l 1199611 . With a 
much bigger sample, we find from SDSS data that ~ 38% 
of all red galaxies in the volume-limited sample have 
[O II] A3727 positively detected. As shown above, the 
red-sequence population actually has two major sub¬ 
components: the LINER/Quiescent sequence and the 
TO/Seyferts. The former dominates the red galaxy pop¬ 
ulation; the LINER-like population is one of the two 
main components in the [O II]-Ha bimodality. Table El 


gives the detailed population breakdown of red galax¬ 
ies in the volume-limited sample. Emission lines are de¬ 
tected in 52.2% of all red galaxies; LINER-like galaxies 
comprise at least 28.8%, dusty-starforming galaxies com¬ 
prise about 6.5%, while TOs and Seyferts are less than 
16.8%. 

4.2. Possible Sourees of the Emission in 
High-[O 11]/Ha Galaxies 

As discussed in the introduction, many recent stud¬ 
ies suggest that star formation could contribute substan¬ 
tially to the [O II] and H/3 emission in Type H QSOs 
and Seyfert 2 galaxies. Is this also the case for LIN- 
ERs? We have shown in ii3.2l that the [O II]-Ha ratio in 
the high-[0 II]/Ha population cannot be reached by star 
formation alone. Here, we provide additional evidence 
using a BPT-like diagram involving [O II] A3727 to show 
that star formation does not contribute substantially to 
[O II] emission in these galaxies. 

Eirst we redivide our emission-line galaxies into four 
categories: star-forming galaxies, Seyferts, LINERs, and 
Transition Objects, using definitions similar to the con¬ 
ventions used in the literature. Our classification scheme 
is illustrated in Fig. M. We firs t emplo y the demarcation 
proposed by iKauffmann et al.1 ll2003ll to separate star¬ 
forming galaxies from other categories. Then we define 
Seyferts conventionally as galaxies with [O III] A5007/H/3 
> 3. The remaining region b elongs to LINERs and TOs. 
We adopt the definition in iHo et a D iIT 9^ to define 
LINERs as objects with [O I] A6300/Ha > 0.17 and TOs 
as objects with [O I] A6300/Ha < 0.17. 

Fig. El shows a BPT-like diagram with 
the [O II] A3727/H/3 ratio plotted against the 
[N II] A6583/Ha ratio. A comparison between panel (a) 
and panel (c) shows clearly that LINER-like galaxies 
have no overlap with star-forming galaxies. If, for 
instance, star formation contributes substantially to the 
[O II] emission in LINERs, it will also contribute to H/3 
in proportion according to a star-forming ratio. As a 
result, the [O II] /H/3 ratio and [N II] /Ha ratio would be 
displaced towards the region occupied by star-forming 
galaxies. In fact, the LINERs have very little vertical 
overlap with star-forming galaxies in [O II] /H/3 ratio. 
Interestingly, Seyferts and Transition Objects do have 
huge overlap with star-forming galaxies, which is con¬ 
sistent with the picture that some of them might have 
substantial fractions of their [O II] and H/3 emission 
contributed by HII regions. Note LINER-like galaxies 
all fall in the high-[0 II] /H/3 category, which is roughly 
equivalent to the high-[0 II] /Ha category. In fact, 
LINER-like galaxies have the highest [O II] /H/3 ratios 
observed. Thus it is impossible to produce LINER-like 
line ratios by combining Seyferts and HII region spectra. 

It is logically possible that pure LINERs have more 
extreme ratios and star formation contribution brought 
them down to the current positions in the line ratio di¬ 
agram. If that were true, the star formation rate would 
have to be remarkably uniform among all LINERs. Oth¬ 
erwise, we would see a spread of emission line ratios, pro¬ 
ducing an elongated distribution like those for Seyferts 
and Transition Objects. Thus this scenario is highly un¬ 
likely. 

It is therefor clear that the high-[0 II]/Ha mode in 
the bimodality is mainly composed of LINERs and 
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Fig. 9.— (a) Line ratio diagnostic diagram (BPT diagram) showing the emission-line flux ratio [O III]/[0 II] versus the ratio [N II]/Ho 
for the magnitude-limited sample (0.05 < z < 0.1). All galaxies plotted have at least 2a detections on all four emission lines. Tradi- 
tional_definitions_forSwf^s and LINERs are indicated by the horizontal and vertical lines |Heckma£j|1980t |Shuder_^^JsterbrocyQ^81J; 

IVeilleux fc Osterbrockl ll987ll . Contours are plotted at number densities of 45,000n galaxies per unit area, where n=l,2,. The effect 

of Ay = 2 of extinction (for Ry = 3.1) is indicated by the arrow, (b) Same diagram as panel (a) but only High-[0 II]/Ha galaxies are 
shown. These are centered in the LINER region. The Icr measurement errors propagated to log-line-ratio is shown as the cross. Contours 
are plotted at number densities of 25,000ri galaxies per unit area, where n=l,2,... . (c) Same diagram but only high-[0 II]/Ha galaxies 
which also have at least 2a detections on [O I] and H/3 are plotted. They have similar distributions as those without detectable [O I] or 
H/3, only biased slightly against high [N II]/Ha ratio and high [O III]/[0 II] ratio. 


Nonzero? 

[O II]/Hq 

Number 

Percentage 

Category 

(O 11] 

Ha 

N 

N 


12913 

47.8% 

Quiescent 

N 

Y 

Uncertain 

3886 

14.4% 

LINERs, TOs and Seyferts ((N II]/Ha > 0.6, 11.1%); 
Dusty SF ([N II]/Ho < 0.6, 3.2%) 

Y 

N 

High 

2209 

8.2% 

Possibly LINER-like 

Y 

Y 

High 

5571 

20.6% 

LINER-like 

Low 

2421 

9.0% 

TOs and Seyferts ([N II]/Ha > 0.6, 5.7%); Dusty SF 
([N II]/Ha < 0.6, 3.3%) 

Total 


27000 

100% 

All red galaxies 


TABLE 2 

Detailed population breakdown of red galaxies in the volume limited sample, see a4.ll 


the line emission is produced by mechanisms other 
than star formation. However, the exact ionization 
mechanism for emissions with LINER-type line ra¬ 
tios is uncertain. The definition of LINER, as in¬ 
dicated by its name, refers only to emission regions 
around the nucleus of a galaxy, but LINER-like line 
ratios have also been observed i n extended em i ssion¬ 
line regions lIPhilliDS et ahl llQSfit iHeckman et alJ 119891: 
[Goiidfrooii et alJ|1994 iZeilinger et alJIlQQbt iGoudfrooiil 
Il999t iSarzi et aljl200,^ Graves et al. in prep). Since 
SDSS spectra were taken with 3” fibers, which cover a 
region ~ 4kpc/h in diameter at the median redshift of 
our sample, both the nuclei and possible extended emis¬ 
sion regions in the galaxies would be included. Thus, we 
cannot tell whether the line emission in these red galax¬ 
ies is associated with compact nuclear sources, extended 


regions, or both. 

However, the occurence of LINER-like line ratios in 
SDSS data is ro ughly consistent with results from the 
Paloinar survey ijEjlioDenko fc Sargent! Il98(il: iHo et all 
ll997albH . which was conducted with a similar aper¬ 
ture (2”x4”) but targeted very nearby galaxies, therefore 
probing much smal ler regions '^100 parsec in diameter. 
Using that dataset. [Ho et al ] lll997bD found that ~ 30% 
of early-type galaxies host a LINER. In the present SDSS 
sample, a similar fraction, namely > 29%, of red galax¬ 
ies have LINER-like emission line ratios. We thus ex¬ 
pect that if emission from extended regions dominates 
the lines we observe, its line ratios and strengths must 
still correlate well with emission from the galaxies’ nuclei. 

Besides photoionization by a central non-stellar source, 
such as an accretion-powered system llFerland fc Netzeil 
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F ig. 10.— This diagr am shows the definition scheme described 
in i4.2l and used in Fi^. We fi rst employ the demarcation pro¬ 
posed by IKanftmann et ahl 1200,^11 to separate star-forming galax¬ 
ies from other categories. Then we defin e Seyferts convention¬ 
ally as galaxies with [O III] A500‘^ H/3 > 3 iShuder Kr. Osterbrocld 
119811; IVeilleux fc Osterbrocki 119871) . Last, we adopt the defini¬ 
tion in Tn^^^^ .1 <1997ali to separate LINERs and TOs using their 
[O I] A6300/Hq; ratio. 


Il98.3t iHalnern fc Steinerl Il983t iGroves et alJ l2004bll . 
many ionization mechanisms have been proposed to 
explain the LINER- and TO-like optical emission in 
early-type galaxies: (1) phot oionization by the X- 
rays radiated as hot gas cools llVoit fc Donahuj 119901: 

iDonahiie fc Voit|ll991llK^ ll989t) , (2) heat tr ansfer from 

hot gas to cooler gas llSoarks et al.111 989H. ISI colli- 
sional excitation by shock waves llHeckman et all 119891: 
iDopita fc Sutherland! Il99,^ . (4) photoionization from 
some comp onent of the stellar population , such as post- 
AGB stars lldi Serego Alighieri et al.l[l99r)HBinette et all 
imi. We have compared the line ratios in the LINER- 
like SDSS galaxies with theoretical predictions o f both 
shockwave models lIDopita fc Sutherland! I1995L I199H1 
and the lowest-ionization Seyfert photoionization mod¬ 
els llGroves et al ] l2nn4albll . Both models produce line 
ratios roughly consistent with the observed ones. De¬ 
tailed comparisons also require a careful extinction cor¬ 
rection (given the possibility of patchiness) and metallic- 
ity measurements, which are beyond the limits of this pa¬ 
per. To reliably distinguish these different mechanisms, 
more information than just line ratio is needed, such as 
the spatial extent and kinematics of the emission-line re¬ 
gions, the spatial variation of line ratios, correlation with 
X-ray, UV and radio properties, etc. It is also quite pos¬ 
sible that many of these mechanisms are involved. This 
is beyond the scope of the current paper, and as such we 
do not try to answer this question here. 

5. IMPLICATIONS FOR STUDIES OF 
POST-STARBURST GALAXIES 


Ha luminosity, the best optical indicator of star 
formation rate, is very difficult to measure at 
z > 0.5 due to atmospheric emission in the near 
infrared. I nstead [OII] A3727 lum i nosity is com¬ 
monl y used lIDressler fc GunnI |1982_|Zabludoffet^y 

1996j__|Hamme^etay__1199?L_|Ii2S£,§L^iil__li^2i 

Dressler^^^ ]1999yPoggianti T_al [1999yBaloehet^y 

1999|__|Rosa^Goi^alezet^ 120021 120031 

Hip pelein et al I2003t lYang et ^ 12004 iTran et alJ 
2004^ . However. we have shown above that the origin of 
[O II] in most red galaxies is not star formation. This 
could have strong consequences for the definition of 
post-starburst galaxy samples, as they are approaching 
t he red sequence and co uld suffer the same issues. 

iQuintero et al.l ll2004ll have demonstrated that post- 
starburst galaxies in the SDSS sample comprise a dis¬ 
tinct population in a plot of Ha EW against A/K ratio'’’. 
The Ha EW provides a measure of the current specific 
SFR in a galaxy, while the A/K ratio is linked to the av¬ 
erage specific SFR over the past 1 Gyr; comparing these 
two quantities indicates how the SFR is changing. In 
the top panel of Fig. d, we reproduce Fig.3 of Quin¬ 
tero et al. with our own measurements for all galaxies 
that are within 0.07 < z < 0.1 and have both [O II] and 
Ha visible in their spectra. For most galaxies, the two 
star formation indicators are highly correlated, as one 
would expect. Post-starburst galaxies stand out as the 
horizontal spur with nearly zero Ha EW but relatively 
large A/K, i.e., they have little or no current SFR but 
large recent SFR. Alert readers will notice the difference 
between our A/K measurement from that of Quintero 
et al.’s; our linear decomposition is done after the sub¬ 
traction of a continuum, while theirs is done without 
continuum subtraction; the templates used also differ. 
In spite of this, the post-starburst population remains 
identifiable. We define post-starbursts as galaxies satis¬ 
fying three criteria: A/K > 0.25, Ha EW < 8(A/K) and 
Ha EW < lOA. 

If we replace Ha by [O II], as shown in the bottom 
panel of Fig. d the post-starburst spur disappears. Al¬ 
though they have minimal Ha emission, post-starburst 
galaxies span a wide range of [O II] EW. In a plot of [O II] 
EW vs. A/K ratio, post-starburst galaxies identified us¬ 
ing Ha overlap with the region occupied by star-forming 
galaxies, as seen in the lower panel of Fig.lT^ 

The left panel of Fig. ^1 shows more directly what 
causes the difference between the two panels of Fig. 1121 
This figure replots Fig.|21but using only galaxies with an 
A/K greater than 0.2 5. This corresponds to the A/K cut 
IQuintero et alJ l|2004ll used to define their post-starburst 
(K-l-A) sample; in this plot, we include all galaxies re¬ 
gardless of Ha EW. As expected, most galaxies on the 
star-forming sequence from Fig.|21remain, while most red 
galaxies are excluded. However, there is a small group of 
points forming a vertical spike near zero Ha EW; these 
are post-starburst galaxies having high A/K and little or 
no Ha emission. There is a gap in Ha that naturally sep- 

A/K ratio is the ratio between the coefficient of the young 
component to the coefficient of the old component used in the linear 
decomposition of a galaxy spectrum, as described in Appendix^ 
a young stellar population will have high A/K ratio 

® In this section, we plot all EW measurements regardless of the 
signal-to-noise of the detection, since galaxies without detectable 
emission are prime candidates to be post-starburst galaxies. 






























































14 Yan et al. 



-1,5 -1.0 -0.5 0.0 -0156 -1.0 -0.5 0.0 -0156 -1.0 -0.5 0.0 0.5 

Log([Nll] A6584/Ha) Log([Nll] A6584/Ha) Log([Nll] X6584/Ha) 


Fig. 11.— BPT-like diagrams showing [O II] A3727/H/3 plotted against [N II] A6583/Ha. The positions of star-forming galaxies (panel 
a), Seyferts (b) and LINERs (c) are shown separately. The effect of Ay = 2 of extinction (for Ry = 3.1) is indicated by the arrow in panel 
a. There is very little vertical overlap between LINERs and star-forming galaxies in (O II]/H/3 ratio, which shows that star formation do not 
contribute substantialy to the (O II] and H/3 emission in LINERs. Note LINERs all fall in the high-(0 II]/H/3 (equivalent to high-(0 II]/Ha) 
category. 


arates this post-starburst population from star-forming 
galaxies. However, there exists no corresponding gap in 
[O II] EW. 

In fact, roughly 50% of post-starburst (K-l-A) galaxies 
identihed using Ha belong to either the high-[0 II]/Ha 
population or the [O II]-only population defined in [}31 
the remainder are classified as quiescent. Although they 
show minimal Ha emission, they have LINER/AGN- 
like line ratios, and can have appreciable [O II] emission. 
Thus, Ha can be used as a universal star formation in¬ 
dicator for most cases with little contamination. How¬ 
ever, in red and post-starburst galaxies, LINER/AGN 
can possess [O II] lines as strong as those in star-forming 
galaxies. 

For this reason, using [O II] to identify post-starbursts 
is problematic. Adopting a low [O II] EW cut would 
make the selection heavily incomplete, while a high cut 
would bring huge contamination from star-forming galax¬ 
ies. In addition, differing levels of measurement errors 
will cause varying degrees of contamination. This could 
well explain the drastically different post-s tarburst abun¬ 
dances found in some previous studies I Drgsslsr^ iiO 
Il999t iPog-gianti et al.l[l999t iBalogh et al.1 19991) . In the 
volume-limited sample, 66% of post-starburst galaxies 
have detectable [OII] emission; 57% have [O II] EW 
greater than 2.5A. Previous studies which used an [O II] 
EW cut of this level could be missing half of all the post- 
starbursts. 

Unlike [O II], Ha appears to provide a generally reli¬ 
able indication of the presence of star formation in the 
absence of a full BPT diagnosis. However, at z > 0.4, Ha 
is difficult to measure due to strong atmospheric emission 
and declining instrumental sensitivity. [O II] has been 


used as a proxy for measuring SFR at redshifts up to 
z ~ 1.4, but as shown above, it is problematic for red and 
post-starburst galaxies. In its place, we suggest the use 
of H/3, which is available out to higher redshifts than Ha, 
but is much more reliable for red gala xies than [O II]. I n 
the DEEP2 Galaxy Redshift Survey llDavis et al ] |200l . 
H/3 can be measured up to z ~ 0.9, enabling a direct com¬ 
parison with low-z studies. Although it is much weaker 
than Ha, more affected by dust, and has a deeper stel¬ 
lar absorption component to remove, H/3 behaves simi¬ 
larly to Ha in SDSS data, as shown in the right panel 
of Fig. ini In a future paper, we will use H/3 to iden¬ 
tify post-starburst galaxies in DEEP2 and compare with 
low-z results. Finally, for completeness, we also give an 
empirical bimodality demarcation for [O II]-H/3: 

EW([0 II]) = I8EW(H/3) - 6. (7) 

5.1. The Nature of Post-starburst emission 

Tradiationally, post-starburst galaxies have been de- 
fined as galaxies w ith no detectable line emission. 
iQiiintero et alJ (j2nn4fl enlarged the definition to include 
galaxies with nonzero Ha emission but too weak com¬ 
pared to star-forming galaxies with a similar fraction of 
young stars (similar A/K ratio, see Fig. 11211 . As shown 
above, these galaxies can have rather strong [O II] lines, 
even when Ha is not detectable (consistent with zero). In 
the volume-limited sample, 78% of post-starburst galax¬ 
ies have at least one emission line detected at > 2a sig¬ 
nificance. What is the nature of this emission? Does it 
indicate residual star formation? 

For post-starbursts that have reliable line detections, 
we can again employ BPT diagrams to study the nature 
of their emission. Fig. m shows two BPT diagrams for 
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Fig. 12.— Top panel: The distribution of Hq EW vs. A/K ratio for all galaxies with 0.05 < 2 : < 0.2 in the SDSS DR4 main galaxy 
sample that have both [O II] and Ha covered in their spectra (compare to Fig.3 in Quintero et ah). Post-starburst galaxies comprise a 
distinct population and stand out as the horizontal spur with nearly zero Ha EW but relatively large A/K ratio. The empirical selection 
criteria are shown with the two solid lines. Bottom panel: The distribution of [O II] EW vs. A/K ratio for the same sample. The horizontal 
spur apparent in the top panel disappears. The dark points indicate the positions in the bottom panel of post-starburst galaxies selected 
in the top panel. These galaxies span a wide range in [O II] EW and overlap with the region occupied by star-forming galaxies. 


galaxies with A/K greater than 0.25 in the magnitude- 
limited sample (0.05 < z < 0.1), i.e., only galaxies which 
contain young stellar populations. As before, we require 
a minimum detection of 2cr on all lines, which excludes 
some fraction of the post-starburst galaxies from the di¬ 
agrams depending on the lines used. For example, about 
1/3 of all post-starbursts are included in the bottom 
panel but only 5% in the top panel. As shown in this 
figure, most post-starburst galaxies in which all emission 
lines involved are positively detected show AGN-type line 
ratios and cover all categories of AGN - LINER, TO, and 
Seyfert. Only 6% of these post-starbursts show line ra¬ 
tios characteristic of star-forming HII regions, which also 
turn out to be preferentially high-Ha EW cases. 

Those post-starbursts that cannot be plotted on the 
BPT diagram mostly belong to the high-[0 II]/Ha group, 
the [O II]-only group, or the quiescent group. Of the 
~ 20% low-[0 Il]/Ha and Ha-only galaxies, the majority 
have [N II]/Ha> 0.6 suggesting a non-starforming origin 
for the emission lines. 

We can break down the classifications for all post- 
starbursts following the scheme used in Table [21 and 
use the [O II]-Ha demarcation and [N II]/Ha ratio to 


discriminate between star-forming galaxies and other 
groups. We find only 5% of post-starbursts in the 
volume limited sample have emission lines dominated 
by residual star formation. Instead, the emission in 
most post-starbursts has AGN-like line ratios. Re¬ 
cently, much evidence has suggest ed a connection 
between AGN and post-starbu rsts |lBrQth£rtQn_£LalJ 
T99g|jGaTia^^ y_ Stqcktonll2nniUKa,iiffmam'^tani2flfl^ 

Hecfanail2004lHc32005ir i)o AGN play a role in quench¬ 

ing star formation? The fact that post-starbursts have 
AGN-like line ratios is consistent with this picture. 

Gonservatively speaking, without knowledge of the 
spatial distribution of emission and information from 
other wavebands, we cannot conclude definitely that 
these are true AGNs. Unlike most red, elliptical galax¬ 
ies, post-starburst galaxies have young stellar popula¬ 
tions, which allow at least one more candidate excitation 
mechanism for LINER- tvpe line ratios aside f rom those 
described above in M4.2I iTaniguchi et al.1 (1200011 point out 
that hot planetary nebula nuclei formed 100-500Myr fol¬ 
lowing a starburst could provide the ionizing spectra nec¬ 
essary to form LINER-type emission lines. Further ob¬ 
servations of post-starburst galaxies in other wavebands 
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Fig. 13. — Left: [O II] EW vs. Ho EW distribution for galaxies with AjK > 0.25; i.e., those with relatively young stellar populations. 
Compare to Fig. Hi. However, here we include all EW measurements regardless of the signal-to-noise of the detection. Right: same galaxies 
as in left panel, but with H/3 EW on the x-axis. In the H/3 distriubtion, although EWs are smaller and measurement errors are larger than 
for Ho;, the separation between post-starburst and star-forming galaxies remain apparent. 


and more detailed modeling should provide a final answer 
to this question. If AGNs are in fact ubiquitous among 
post-starburst galaxies, it may provide a hint that AGN 
played a role in the quenching of their star formation. 

6. SUMMARY 

We have measured the fluxes and EWs of [O II], Ha 
and several other emission lines for ~ 300, 000 galaxies 
in the SDSS DR4 main galaxy sample after careful sub¬ 
traction of the stellar continua. From the comparison 
between [O II] and Ha, combined with a variety of line 
ratio diagnostics, we have investigated the origin of line 
emission in red galaxies and the implications for post- 
starburst galaxy studies. Our main conclusions are as 
follows: 

1. Galaxies display a bimodality in [O II]/Ha ratio, 
which corresponds closely to the bimodality in their 
rest-frame colors. In blue galaxies, the [O II]/Ha 
ratio generally matches the expectation for star¬ 
forming HII regions. However, in most red galaxies, 
the JO II]/Ha ratio is unusually high and difficult 
to reconcile with predictions for star-forming HII 
regions. 

2. About 52% of all red galaxies have detectable line 
emission; 38% have detectable [O II] A3727 emis¬ 
sion. More than 29% of all red galaxies have line ra¬ 
tios characteristic of LINERs, while less than 17% 
are TOs or Seyferts, and ~ 6% have lines domi¬ 
nated by star formation (probably dusty starbursts 
and edge-on spirals). Further information about 
the spatial distribution of emission in these galaxies 


and/or X-ray, UV and radio observation is neces¬ 
sary to firmly identify the origin of the line emission 
in each case. 

3. LINER-like galaxies make up the high-[0 II]/Ha 
mode in the [OII]/Ha bimodality. Other 
than modest differences in the luminosity range 
spanned, they are essentially indistinguishable 
from quiescent galaxies in the color-magnitude- 
concentration space. The combination of LINER- 
like and quiescent galaxies defines a uniform red 
sequence in color-magnitude-concentration space, 
which can be effectively selected by a simple divi¬ 
sion in the [O II]-Ha EW diagram. The remain¬ 
ing red galaxies, those with low [O II]/Ha ratio, 
are mostly TOs, dusty-starforming galaxies, and 
a small fraction of Seyferts. This group also have 
lower galactic concentration and slightly bluer color 
than the LINER/Quiescent red sequence. 

4. Post-starburst galaxies, identified in the SDSS 
dataset using a lack of Ha emission to indicate that 
star formation has ceased, often exhibit significant 
]0 II] EW. Their positions in the [O II]-Ha plot 
and BPT diagrams are the same as for Seyferts, 
LINERs and TOs, which suggest they may harbor 
AGNs. Less than 5% of this post-starburst sample 
show evidence for residual star formation. 

5. (O II] emission in red galaxies and in post-starburst 
galaxies (i.e., in cases where it is not induced by 
star-formation) can be as strong as in star-forming 
galaxies. Locally at least, [O II] can only be used 
as a star-formation indicator for blue galaxies. 
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recommend using H/3 as an alternative for high- 
redshift studies. 
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APPENDIX 

STELLAR CONTINUUM SUBTRACTION 

To get an accurate measurement of the emission line fluxes and equivalent widths, we need to subtract the stellar 
component underneath. All Balmer lines have a broad absorption component, while the continuum underlying [O II] 
has many narrower absorption features. The best way of removing the stellar component is often to fit the entire 
spectrum with a stellar population synthesis model. Here, we trade off this freedom for robustness and efficiency. 

Instead of doing full population synthesis modeling, we decompose the spectra into two components: an old stellar 
populatio n and a young s t ellar p opulation. Similar methods have been used in modeling the spectra of E-l-A galaxies, 
such as in iQiiintero et al.1 ll2004ll. who employed an A -star spectrum and a composite luminous red galaxy spectrum as 
templates. As pointed out bv iDressler V Gunul ill 98,111 . a linear combination of the two components can simultaneously 
reproduce both the broad spectral continuum and the depths of narrow features for post-starburst galaxies, over the 
range of 3400 to 5400 angstroms. For galaxies other than post-starbursts, the broadband continuum is not always 
matchable by this method, but the depths of the narrow features can always be well fitted by a combination of these two 
components; this is sufficient for our purposes. Thus, we model the continuum-subtracted spectra by the combination 
of the two templates, each continuum-subtracted in the same way. The broadband continuum subtracted is determined 
by the inverse variance-weighted mean of the spectrum in a 300A window. The method works well for the regions 
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Fig. 14.— Line ratio diagnostic diagrams for galaxies with 
AjK > 0.25 in the magnitude limited sample: [O III] A5007/H/3 
vs. [N II] A6583 /Hq: (top) and [O III] A5007/[0 II] A3727 vs. 
[N II] A6583 /Hq; (bottom). All galaxies plotted here have at least 
2cr detections on the emission lines involved. Post-starburst galax¬ 
ies are highlighted with the diamond symbol; they mostly have 
LINER-type line ratios. The number of post-starburst galaxies 
plotted in the top panel is much smaller because H/3 is more diffi¬ 
cult to detect than the strong lines used in the bottom panel. 


6. More than half of all the post-starbursts have de¬ 
tectable [O II] emission. Post-starburst samples 
defined using [O II] as a star-formation indicator 
will be very incomplete, especially if AGN/LINER 
rates and intensity were higher in the past. We 
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around all emission lines used here: viz., [O II] A3727, H/3, [O III] A5007, JO I], Ha, and [N II] A6583; see iiA.2l 

The coefficients of the two templates are obtained from a linear fit over the restframe wavelength range 3600A < 
A < 5400A. Before fitting, both templates are convolved with a Gaussian kernel with a dispersion equivalent to the 
combination of the instrumental dispersion and the velocity dispersion of the galaxy being fitted in quadrature. The 
fitting is weighted by the inverse variance in the observed spectrum. We mask regions around emission lines using 
adaptive windows set by the velocity dispersion measurement of each galaxy. 

Spectra that cover a total of less than 900A (3 times the broadband continuum subtraction window) in the observed 
frame or less than 200A out of the restframe wavelength range from 3600A to 5400A are excluded from the samples 
used in this paper, in order to ensure good removal of stellar continuum features. 

The Templates 

Our young stellar popul ation template is given by a model spectrum built using the new population synthesis code of 
iRruzual Kr. Cha.rk^ 112(10311 . The template spectrum is taken at 0.3 Gyr after a starburst with a constant star formation 
rate lasting 0.1 Gyr with solar metallicity. This template is used instead of an A-star spectrum because it is a more 
physical choice. Moreover, it includes all the spectral features produced in intermediate and low mass stars which also 
exist in a starburst. This particular choice of spectrum also provides a relatively good subtraction of the continuum 
features for galaxies with ongoing star formation. 

The old stellar population template is also built using iBruzual fc Gharlod ll2003ll code. The model used is a 7Gyr 
old simple stellar population with solar metallicity. On e could have instead u sed a purely empirical template, such 
as the SDSS Luminous Red Galaxy coadded spectrum ijEisenstein et a,l1l2flfl3fl . This spectrum closely resembles our 
synthesized old stellar population spectrum if the broadband continuum is subtracted: almost every feature in the two 
spectra matches. However, such a spectrum appears to include the sort of line emission we are looking for; e.g., it has 
an [O II] EW of ~ 2A. Thus, we avoided using it as a template. 

Testing the subtraction 

We have tested our subtraction of the stellar continuum using a set of population synthesis models and real spectra. 
Except for models with very recent bursts (younger than 200 Myr), the residuals of the fits are sufficiently small (the 
residual EWs in H/3 are within ±0.4A and those in Ha are within ±0.7A). For recent bursts the discrepancy occurs 
because such galaxies are likely to be dominated by B stars, which have narrower Balmer absorption lines than A 
dwarfs. A characteristic residual is therefore apparent for these galaxies. However, we have found this sort of feature 
only very rarely (< 1%) in visual inspection of residuals from a subset of the spectra used here. Since this paper 
focuses on emission lines in red galaxies, the poor residuals for extremely young bursts do not affect our conclusions. 
In some cases, there could also be systematic residuals around the Ca K line, which is a metallicity-sensitive feature. 
Since this f eature does not overlap any of the emission lines we measure, it does not affect our results. 

Fig. lAI 5l shows three examples of the stellar continuum subtraction. They are selected to have a S/N ratio near the 
median value in the sample. The residual of the subtraction is shown at the bottom of each panel. Most stellar features 
are successfully subtracted with small residuals. In high S/N spectra where residuals are not dominated by noise, the 
median root-mean-square of the residual around each line is less than 2% of the continuum level except around [O II], 
where residuals approach 10% (this is an uppper limit, as noise in the spectra is higher at shorter wavelengths). 

LINE MEASUREMENTS 

We measure the emission line flux by two methods: flux summing and Gaussian fitting. Table IB41 lists our line 
definitions. For each emission line, we define a center passband in which we measure the line flux, and two sidebands 
in which we measure the continuum flux. We use these same wavelength ranges to establish which spectra are well- 
measured in the vicinity of Ha and [O II], defining our sample. We require that at least 83% of all pixels in the center 
passband and 17% of the pixels in each sideband be well-measured for both [O II]and Ha. In order to be considered 
good, we require both a pixel and its four adjoining pixels to have inverse variance above the greater of and 

one-ninth the median inverse variance for the spectrum. 

First, the stellar continuum is reconstructed via our linear fit with two templates. Although the fit is performed with 
the broadband continuum subtracted templates, in reconstruction, we use the original templates with the two linear 
coefficients from the fit. This reconstructed stellar continuum has correct depths for narrow features but offset overall 
amplitude. We then add a constant to it to match the median of the observed flux in the two sidebands, respectively 
for each emission line considered. Subtracting off this stellar continuum from the observed spectrum then leaves us 
with an emission-line only spectrum. The line flux is measured by integration in the center passband. 

In addition to measuring total flux from simple flux-summing in a window, we also apply a Gaussian fitting technique. 
Here H/3, [O III] A5007, and [O I] A6300 are each modeled as single Gaussians, where the line amplitude, dispersion and 
central wavelength are fit simultaneously by a Levenberg-Marquardt least-squares method with realistic boundaries 
set on the last two parameters. [O II] A3727 is fitted as a doublet with the wavelength ratio fixed and dispersion shared 
between the two lines. Because of their proximity, Ha and the [N II] AA6548,6583 doublet are fit simultaneously with 
three Gaussians with wavelength ratios fixed and dispersion shared. We have checked the line fluxes obtained from 
flux-summing against those from Gaussian fitting. Good agreement is achieved for lines with detectable EWs. For 
instance, the standard deviation of the fractional difference in Ha flux is about 7% for galaxies with Ha EW around 
lA; that of [O II] flux is about 16% at EW([0 II])~ 4A. Systematic differences occur near zero flux, as the Gaussian 
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Fig. A15.— Sample spectra illustrating the accuracy of our stellar continuum subtraction. These spectra are selected to have a S/N 
ratio near the median value in the volume-limited sample. The fit and residual are offset for clarity. Weak emission lines are more easily 
seen in the residual than in the original spectra, such as [O III] A5007 in the top panel and [O I] A6300 in all three spectra. Also worth 
noting is that the different line ratios show the different origins of the line emission, as discussed in m and HD 


fitting routine used tends to avoid low amplitudes, fit a broad Gaussian to the noise and overestimate the flux. For 
this reason, we use the flux-summing measurements for all results presented in this paper. 

As described in a and illustrated in Fig. n the measured EW distributions of galaxies that fall to the left of the 
demarcation in the [O II]-Ha EW diagram can be modeled well by the combination of a Gaussian distribution centered 
near zero and a log-normal distribution starting from the peak of the normal distribution; these two components 
correspond very closely to the quiescent population and the LINER-like population, respectively. We make use of 
this fact to measure zero point offsets for our em ission line measurements, as presumably the true EW for quiescent 
galaxies is zero or extremely close to it. Table | bH 1 lists the zero point determined by this method for each emission line 
used in this paper and the a of the Gaussian distribution of that line’s EW for quiescent galaxies. The latter gives a 
measure of the line detectability. We have subtracted the zero point derived by this method from all the EW values 
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Line 

Zeropoint (A) cr (A) 

JO II] A3727 

1.45 

1.56 

H/3 

0.318 

0.231 

JO III] A5007 

0.54 

0.35 

[O I] A6300 

0.119 

0.238 

Hq 

0.487 

0.349 

IN II] A6583 

0.195 

0.312 

TABLE B3 



For each emission line used in this paper, we list the zero point offset of the original EW measurements and the standard 
DEVIATION ((t) in EW FOR GALAXIES WHOSE EW IS CONSISTENT WITH ZERO. THESE VALUES HAVE BEEN DETERMINED FROM THE GAUSSIAN 
+ LOG-NORMAL DECOMPOSITION OF THE EW DISTRIBUTION FOR GALAXIES TO THE LEFT OF THE DEMARCATION IN [O II]-HaEW DIAGRAM. 


Index Bandpass 

[O II] A3727 3716.3 - 3738.3 

H/3 4857.45-4867.05 “ 

[O III] A5007 4998.2-5018.2 

[O I] A6300 6292.0-6312.0 

Ho 6554.6 - 6574.6 

[N II] A6583 6575.3 - 6595.3 


Blue Sideband 


Red Sideband 


Reference 


3696.3 -3716.3 3738.3 - 3758.3 Fisher et al. 11998') 

4798.875-4838.875 4885.625 - 4925.625 ^isheret^ (199^ 
4978.2-4998.2 5018.2-5038.2 

6272.0-6292.0 6312.0-6332.0 

6483.0 - 6513.0 6623.0 - 6653.0 

6483.0 - 6513.0 6623.0 - 6653.0 


‘‘The central bandpass is modified to cover only the middle 1/3 of the window in the|Fisher_et_alJ (1993) definition. The integrated flux in 
the modified window agrees better (i.e., has less scatter) with the Gaussian fitting results than the flux in the wider, original window does. 


TABLE B4 

Window definitions for the emission lines we measure. 


used in this paper. The cr’s of the Gaussian distributions for Ha and H/3 suggest the EW errors for these two lines are 
too low by 33%, perhaps due to the difficulty of subtracting broad Balmer absorption features. We thus multiply the 
errors used for these lines in the paper by 1.5 when determining significances, etc. For all other lines, the Gaussian 
width is in good agreement with the estimated errors. 

We have also checked our EW measurements against those measured by Ghristy Tremonti Isee lTremouti et alJ2f)n4j) . 
In almost every case, any zeropoint offset is small enough to be ignored: the difference be tween the two measure ments 
are of order of the error estimates. However, for H/3, our EW is about lA smaller than iTremonti et al.l ll2004l) when 
H/3 is greater than lOA. This is of little concern, because our main interest in this paper centers around red galaxies 
with little H/3 emission, where the zeropoint offset is less than O.lA. The difference is likely to be due to differences 
in templates and methods for stellar continuum subtraction; galaxies with strong star formation will typically exhibit 
both strong H/3 emission and significant stellar H/3 absorption features. 
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